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Abstract 

We review some results obtained by us and others concerning the structure 
of higher order perturbation series in perturbative QCD and their resummation 
using the renormalization group equation. We illustrate our results by a number of 
examples involving hadronic r decays and e^e~ annihilation. 



* Contribution to the Festschrift published on the occasion of the 60th birthday of Dmitri Kazakov 



Contents 

1 Introduction 3 

2 High order perturbation theory in QCD 5 

2.1 Comparison of observables in e"^e~ annihilation and r decays 5 

2.2 Signal of asymptotic behaviour in r moments l7 

2.3 Us from the r width in a RG invariant way llO 

3 Resummation on the cut: > |l2 

3.1 Example with an explicit solution for q;s(/U^) [12 

3.2 Example with an explicit expression 

for the integral over the effective charge [13 

4 Resummation on the contour |l4 

4.1 Resummation using the four- loop /3 function 15 

4.2 Relations between observables using resummation [l5 

5 Comparison of resummation techniques [16 

6 Summary |23 
A Appendix [24 



2 



1 Introduction 



The mathematical framework of Quantum Field Theory (QFT) is a powerful tool for the 
description of elementary particle interactions. The concept of renormalization and the 
technique of the -R-operation generate renormalizable models of QFT as self-consistent 
tools for a phenomenological analysis [1]. The Standard Model (SM) of "all interactions" 
is formulated entirely as a local quantum field theoretical model. The proof of its renor- 
malizability has required some extension of the standard techniques [21 13]. Despite the 
success in explaining the properties of interactions at high energies, a closed-form solution 
of the SM is still missing. The main route for a quantitative investigation is still the use 
of perturbation theory (PT) where one expands in the small coupling constant where the 
coefficients of the perturbation series are computed in terms of integrals represented by 
Feynman diagrams. 

The classical example of the perturbation theory approach has been developed in the 
context of QED which is part of the SM. After renormalization of QED there are two 
parameters - the electron mass and the fine structure constant a - which are to be fitted 
to observables in terms of a perturbation series in a. With a = 1/137 being small, the 
precision of the results is impressive. However, the question has been raised about the 
behaviour of PT series at high orders. Arguments have been put forward that the series 
is merely an asymptotic one |H |5]. This has caused much interest in the summation 
of asymptotic series in QFT [6]. Different theoretical aspects of the problem have been 
discussed in Ref. [7]. A quantitative investigation of the perturbation series in QFT and 
the question of its convergence has been presented in Ref. [8]. 

As concerns the problem of PT convergence it has been realized that one can identify 
some "important" terms in the entire PT series and resum them. For example, "large logs" 
have first been resummed in [9]. These first results even provoked extreme statements 
about the fate of local QFT [10] . The correct interpretation of the "resummed" results 
emerged within the technique of the renormalization group [TT] (as a review, see Ref. |12j). 
Different types of resummation have been suggested later on: the analysis of classical 
solutions [13], the famous example of the 1/Nc expansion in QCD [13], the large Pq limit 
in QED [15], and the naive nonabelianization in QCD |16j . 
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A precision analysis of the SM requires an extensive use of the renormahzation group 
(RG). The energy scales of SM are quite different and the coupling constants (especially 
as) are quite large. Indeed, the electroweak scale given by Mz,w ~ 90 GeV and nit = 
175 GeV is quite high while the hadronic scales are rather low as e.g. given by rrih = 5 GeV 
and the light hadron masses at around 1 GeV. The process of deep inelastic scattering 
allows one to probe a wide range of scales by varying Q. Modern experiments include 
hadrons (LHC) that requires the analysis of hadronization and the use of QCD at this 
scale with ^^(l GeV) = 0.45. Thus the use of PT in ^^(/i) at low energies /i = 1 GeV 
requires resummation as the convergence is slow. 

Let us list a few technical issues used in the perturbative description of the phe- 
nomenology of the SM: 

i) Use of higher order PT corrections including NNLO corrections which has almost 
become a standard; 

ii) Precise definition of the expansion parameter or the choice of the renormahzation 
scheme. The MS-scheme is a standard for technical reasons since dimensional reg- 
ularization has become the main computational framework for the calculation of 
radiative corrections. However, other more physical schemes are also in use such as 
ay from the Coloumb part of the potential of heavy quarks when using nonrela- 
tivistic QCD to describe the dynamics of heavy quark production near threshold; 

iii) Resummation of some infinite subsets of the perturbation series in a^. 

The most popular way of improving perturbation theory is the use of the RG to sum 
powers of e.g. logs (as ln((5//i))", as(asln((5/yu))". But some other ways are also used: 
the Po dominance (naive nonabelianization) which sums terms of the form {Poas)"', or 
effects of analytic continuation between Euclidean and Minkowskian regions that basically 
deal with terms (vr/Joas)^" [HI CHI [19] . 

In case of using infinite subsets of the perturbation series, PT allows for nonpolynomial 
terms in a^. The PT series are asymptotic and resummation may provide terms that 
interfere with non-perturbative expansions for the description of some processes that 
account for terms of the form exp —{l/as{Q))'- higher twists in light-cone type expansions. 
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or condensate type terms for the expansion of correlation functions at short distance. This 
means that numerical values of such non-perturbative parameters depends on how the PT 
series are treated. 

This is very important for the analysis of hadronic r decays. The theoretical descrip- 
tion is simple and related to e~^e~ annihilation. An unprecedented number of PT terms 
are available. It is clear that the structure of the series is important. Also it is not an 
academic exercise as it is important for hadronic contributions to aEM{Mz) [201 El] 
to the muon g — 2 [22| 123] which are important for constraints on the Higgs mass and are 
also key players in constraining new physics search beyond the SM [2l]. And, it is also a 
field of interest of Dima Kazakov. 

2 High order perturbation theory in QCD 

A classic example of application of the PT series in QCD is the analysis of the e~^e~ 
cross section [251 |26l |27]. The relevant expansion parameter is large, such that a requisite 
accuracy requires a high order PT expansion. With such a large number of terms one 
may already encounter the asymptotic nature of the perturbation series in which case no 
further increase of precision is possible. The main problem for the theory is convergence 
and the interpretation of the numerical values given by the series. An additional freedom 
and also complication is that the expansion parameter is not uniquely determined and the 
series should be analyzed in a scheme- invariant way [28] [29] . Because of the freedom to 
redefine the scheme it is difficult to judge the quality of convergence of the series. In this 
section we present a way to bypass this complication by establishing a relation between 
observables. 

2.1 Comparison of observables in e+e annihilation and r decays 

Within massless perturbative QCD the same Green's function determines the hadronic 
contribution to the r-decay width and the moments of the e~^e~ cross section. This 
allows one to obtain relations between physical observables in the two processes up to an 
unprecedented high order of perturbative QCD [30] . A precision measurement of the r 
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decay width allows one then to predict with high accuracy the first few moments of the 
spectral density in e~^e~ annihilation integrated up to s = m^. 

The question of numerical convergence is influenced to a large extent by the freedom 
of choosing the renormalization scheme for the truncated perturbation series [3T | [32 | [33]. 
Therefore, it is desirable to obtain predictions for observables which are renormalization 
scheme independent. 

We compare moments of the spectral density in e~^e~ annihilation and the hadronic 
contributions to r(r — )■ u^- + hadrons) [SU [351 IMl EZ]. The reduced decay width 
appearing in 

i^.= ^,^;-^'^- + ^"^^""^^ =3(l + r.) (1) 

is determined by massless perturbative QCD, for which the axial and vector contributions 
are identical. The expansion for r^- starts directly with alfi"^), where a = agju. In 
e^e~ annihilation the cross section is determined by the imaginary part of the vacuum 
polarization function, 

= 127rJmn(s) = iVeJ]Q?(l + r(s)) = 2(l + r(s)). (2) 
In perturbative QCD one has 

r(s) = a(/i2) + (A;i + /3oL)a2(/i2)+ (A;2-^7rX + (2/3oA;i + /3i)L + /3o'^')a'(yt/'^ 
+ (A;3 - 7r2/3o'A;i - ^-ix'^Mi + (3/3oA;2 + 2/3iA;i + - tt^/^^)/, 
+ (3/3oA;i + \fii)L^ + filL%\^^) + . . . (3) 



2 

with L = ln(/i^/s). We define moments of r{s 



rn{so) = {n+l) / — - r{s) (4) 



■^o V'^o 
such that 

r^- = 2ro{ml) - 2r2{ml) + rs{ml). (5) 

Eq. ([5]) can be inverted within perturbation theory. One can then express the perturbative 
representation of one observable, i.e. any given e~^e~ moment r„(m^), in powers of r-r using 
the perturbative expansion of the r-decay observable. The strong coupling constant 
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in any given scheme serves only as an intermediate agent to obtain relations between 
physical observables. The reexpression of one perturbative observable through another is 
a perfectly legitimate procedure in perturbation theory and the result is independent of 
the choice of the renormalization scheme. One finds 

rn{ml) = UnTr + fmr^ + f2nrl + fsnrt + Un4 + 0(rf), (6) 

where the coefficients are given in the Appendix. 

For r°^P = 0.216 ± 0.005 one can investigate the convergence properties of the series 
for the first few moments. One obtains 

ro/0.216 = 1 - 0.284 - 0.069 + 0. 110 + ... (7) 

ri/0.216 = 1 - 0.527 - 0.143 + 0. 177 +... (8) 

rs = 0.216(1 - 0.608 - 0.115 + 0.269 + ... (9) 

rg = 0.216(1 - 0.648 - 0.091 + 0.317+... (10) 

One clearly observes the divergent behaviour of the perturbation series for the moments 
n = 0, 1, 2, 3. Since the relations between different sets of observables are scheme indepen- 
dent, there is no freedom in the redefinition of the expansion parameter. One sees very 
poor convergence. Therefore, the question is whether the asymptotic regime has already 
been reached. 



2.2 Signal of asymptotic behaviour in r moments 

How can the asymptotic nature of the series reveal itself? With a redefinition of the 
charge one can create any type of the series that hides the true rate of convergence. 
Therefore, one should work in a scheme invariant way. We concentrate on the analysis of 
the r-system only. 

The spectral density has been calculated with a very high degree of accuracy within 
perturbation theory (see e.g. [381 EHl HOI E]) and has been confronted with experimental 
data to a very high degree of precision [Ml |35l [361 EZ] • In Ref . [12] arguments have been 
given that, within the finite order perturbation theory, the ultimate theoretical precision 
has been reached already now. The limit of precision exists due to the asymptotic nature 
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of the perturbation theory series. The actual magnitude of this hmiting precision depends 

on the numerical value of the coupling constant which is the expansion parameter. We 

do not discuss power corrections here |l3llll]. 

The central quantity of interest in the r system is the hadronic spectral density which 

can be measured in the finite energy interval (0, M,- = 1.777 GeV). The appropriate 

quantities to be analyzed are the moments. We define moments of the spectral density 

by ^ 

Mn = {n + l) f p(s)s"ds = l + m„ (11) 
Jo 

{Mt- is chosen to be the unity of mass). The invariant content of the investigation of 
the spectrum, i.e. independent of any definition of the charge, would be the simultaneous 
analysis of all the moments. 

In order to get rid of artificial scheme-dependent constants in the perturbation theory 
expressions for the moments we define an effective coupling a{s) directly on the physical 
cut by the relation 

p{s) = l + a{s). (12) 

All the constants that may appear due to a particular choice of the renormalization scheme 
are absorbed into the definition of the effective charge (see e.g. [311 HSl HHl HZ])- When 
defining the effective charge directly by p{s) itself we obtain perturbative corrections to 
the moments only because of the running of a^. Without ag running one would have 

Mn = l + a{Mr) = l + a or mn = a (13) 

At any given order of PT the running of the coupling a{s) contains only logarithms of s, 

a(s) = a + PoLa^ + (/3iL + f3lL')a^ + (/^si^ + ^A/^o^' + PlL^)a^ + . . . (14) 

where a = a{M^), L = \n{M^/s). At fixed order of PT the effects of running die out for 
large n moments improving the convergence of the series 

mo = a + 2.25a2 + 14.13a^ + 87.66a^ + 654.16a^ 
mi = a + 1.125a2 + 4.531a^ + 6.949a^- 64.77a^ 
m2 = a + 0.75a^ + 2.458a^ - 1.032a^ - 68.98a^ (15) 
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For large n the moments behave better because the infrared region of integration is sup- 
pressed but in high orders they start to diverge. For large n, the coefficients of the series 
in Eq. ( |T5l) are saturated by the lowest power of the logarithm. 

To suppress experimental errors from the high energy end of the spectrum, the modified 
system of moments 

Mki = -^^^ p{s){l - s)''ds = l+mk (16) 

can be used. The integral in Eq. (fT6|) is dominated by contributions from the low scale 
region. A disadvantage of choosing such moments is that the factor (1 — s)^ enhances the 
infrared region strongly and ruins the perturbation theory convergence. As an example 
one has 

mo = a + 2.25a2 + 14.13a^ + 87.66a^ + 654.2a^ 

mi = a + 3.3750^ + 23.72a^ + 168.4a^ + 1373.29a^ (17) 

which shows a poor convergence. The reason is the contribution of the logarithmic terms 



and 




which grow as ln(A;) and ln^(A;) for large k. 

The large difference in accuracy between the moments mo and mi is a general feature 
of the moment observables at the fifth order of perturbation theory: one cannot get a 
uniform smallness at this order for several moments at the same time. For any single 
moment, one can always redefine the charge and make the series converge well at any 
desired rate, but then other moments become bad in terms of this redefined charge. The 
invariant statement about the asymptotic growth is that the system of moments m„ 
including m„=o cannot be treated perturbatively at the fifth order of perturbation theory. 

In order to demonstrate this in a scheme invariant way, we choose the second moment 
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(which is aheady well convergent) as a definition of our experimental charge and find 

mo = m2 + 1.5m^ + 9.417m^ + 59.28m^ + 457.54m^, 
mi = ma + 0.375m^ + 1.51ml + 2.527m^ - 17.64m^, 



There is no convergence. 

2.3 as from the r width in a RG invariant way 

Having in mind that the series expansion has reached the ultimate accuracy, we try to 
avoid expansions and analyze the system in a different but still concise way [ISj. The ob- 
servation is that any perturbation theory observable generates a scale due to dimensional 
transmutation and this is its internal scale. It is natural for a numerical analysis (and 
this is also our suggestion) to determine this scale first and then to transform the result 
into a MS-scheme using renormalization group invariance. 

We use the explicit renormalization scheme invariance of the theory to bring the result 
of the perturbation theory calculation into a special scheme first, followed by a numerical 
analysis in this particular scheme. Only after that we transform the obtained numbers 
into the MS reference-scheme. 

A dimensional scale in QCD emerges as a boundary value parameterizing the evolution 
trajectory of the coupling constant. The renormalization group equation 



1712 = 



m2, 



1713 = 



m2 - 0.19m^ - 0.544m^ + 0.742m^ + 16.8m^, 
7712 - 0.3m^ - 0.803m^ + 1.69m^ + 27.2m^. 



(20) 




a 

a = — 
n 



(21) 



is solved by the integral 




) 



<|.(a(/i2)) + 




(22) 



with 




(23) 
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/33 ^ /33 + 4fi:?/3o + 2/t3/3o + fi:?/3i -2fi:i(3fi:2/3o + /32 



The MS-scheme parameter A is defined through the expansion 

The evolution trajectory of the couphng constant is parametrized by the scale parameter 
A and the coefficients of the /3 function /3j with i > 2 (see e.g. Ref. [32]). The evolution 
is invariant under the renormalization group transformation 

a — )■ a(l + Kifl + K2a^ + K^a"^ + . . .) (25) 

with the simultaneous change 

A^ ^ A^e"'^^/*, (26) 

where /3o,i is left invariant and 

This invariance is violated at higher orders of the coupling constant because one has 
omitted higher orders for the /3 functions. This is the source for different numerical 
outputs of analyses in different schemes. 

We introduce an effective charge = [311 ST] and extract the parameter A,- which 
is associated with through Eq. fl22|) with an effective (3 function 

(3^(a^) = -al{2.25 + 4a^ - 12.30^ + 38.1ai^). (27) 

In this procedure the only perturbative objects present are the /3 functions related to 
e^e~ annihilation (in Eq. ( 12T]) ) and r decay (see Eq. ( 127|) ). We treat them both as 
concise expressions, and at every order of the analysis we use the whole information of 
the perturbation theory calculation. For the coupling constant in the MS scheme we 
finally find 

a, = 0.3184 ±0.0159. (28) 

For the reference value of the coupling constant at the scale Mz = 91.187 GeV we run to 
this reference scale with the four-loop (3 function in the MS scheme [19] and three-loop 
matching conditions at the heavy quark (charm and bottom) thresholds [50] to get [H] 

as{Mz) = 0.1184 ± 0.0007exp ± 0.0006^6 (29) 
Can one do better than finite order PT? Yes, but then one has to resum! 
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3 Resummation on the cut: > 



Since new higher order terms in the perturbation series will not become available in the 
near future for e~^e~ annihilation or the r decays, it is tempting to speculate on the 
general structure of the series within PT [511 ESI [S3] • Much attention has been recently 
paid to possible factorial divergences in the PT series [511 EHl EH E71 ES] which is generated 
through the integration over an infrared region in momentum space |59j . 

Going beyond finite order PT directly by taking into account RG logs in p{s) does not 
work. Using a simple approximation p{t) = as{t) and 



with factorial growth which is not Borel summable. The reason is the Landau pole in 
the expression for as(t) or the divergence of the integrand outside the convergence circle 
|a(s) ln(s/t)| < 1. Higher order terms in p{t) are important [60] . 

For an estimate of the uncertainties of the theoretical predictions for the r lepton 
width, different approaches have been used for the definition of the integration over the in- 
frared region. This problem has been widely discussed in the literature (see, e.g. Ref. |61]). 
We propose a set of schemes that regularize the infrared behaviour of the coupling con- 
stant in general and allow one to use any reference scheme for the high energy domain |60] . 
All these schemes are perturbatively equivalent at high energies. The uncertainties that 
come from the low energy region are quite essential as our study will show. 

3.1 Example with an explicit solution for as{fJp') 

Consider first an example with an explicit solution for as{pi^) ^0]. Consider a P function 



with /3"'*(a) = — + . . . at a — 0. The RG equation has a solution 




one finds 
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and the pole at /i^ = of the asymptotic solution a"''(/i^) = (ln(/x^/A^))~^ disappears. 

The example shows that a particular way of summing an infinite number of specific 
perturbative terms for the (3 function can cure the Landau pole problem. No nonper- 
turbative terms are added but the freedom of choosing a renormalization scheme for an 
infinite series was used. This result can be considered either as a pure PT result in some 
particular RGE after an infinite resummation or as a sort of Pade approximation of some 
real /3 function that may also include nonperturbative terms. The only important point is 
that the running coupling obeying the RG equation with such a /3 function has a smooth 
continuation to the infrared region. Because the expansion parameter becomes large in 
the infrared region, the polynomial approximation is invalid in this domain. This is a 
particular example of the general situation that the expansion in the unphysical parame- 
ter as is incorrect and the appropriate way to proceed is to expand one physical quantity 
through another. 

3.2 Example with an explicit expression 

for the integral over the effective charge 

Consider an example with 

= TTYa 

In this case the integral 

1 r 

F(s) = - / a(t)dt 
s Jo 

can be found explicitly [60]. The RG equation for the effective charge a{s) is given by 

ln(s/A2) = -21na(s), 
ays) 

and F{s) reads 

If 1 

F(s) = - / a(t)dt = a(s) + a(sY — a(sY exp( -— ). 

s Jq a{s) 

The last term gives the "condensate" contribution up to logarithmic corrections ~ 
A^/s. This example shows that a change of the evolution in the infrared region resums 
factorials. The last term cannot be detected by integrating the series near a = 0. 
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In order to study higher orders of PT, the Borel transformation is often used. The 
Borel analysis for this example reads 

F{a) = / e~^^''B{^)d^ (30) 
Jo 

and the Borel image B{x) is B{x) = 1 + x + (1 — x)9{x — 1). The PT series with 
all coefficients known does not allow one to restore the exact answer through Borel 
summation. A naive Borel image for polynomials from the definition f l30|) behaves as 
B{x) = Ylk=i fkX^ /{^ ~ 1)' which is correct at small x but not at large x. The explicit 
result shows that the Borel image is singular at x = 1. 

There is still another possibility to treat the PT series: resummation in the complex 
plane [19]. 



4 Resummation on the contour 

By integrating the function T]^^'^{z) over a contour in the complex plane above the 
physical cut s > one finds 

i) 'n.{z)dz = / p{s)ds 

Jc Jcut 

with 

p(s) = -i,(n(s + zO)-n(s-zO)). 

2m 

One then uses the approximation n'^'^*^(2;)|2gc' ~ n^"'"(-2)|2ec which is well justified when 
one is sufficiently far from the physical cut. One then obtains 

had/ \ T„ / „/-„\j„ I ttPT 



Il'''^{z)dz= / p{s)ds= f Il'^'{z)dz 

C Jcut Jc 

where the integral over the hadronic spectrum is computable in perturbative QCD. The 
total decay rate of the r and its moments are the quantities that can be computed in this 
way. The use of the RG improved 11^^(2) on the contour has been first considered in [12]. 
The technique is now known as Contour Improved Perturbation Theory (CIPT) [T9l [62] . 
Parametrizing the contour by = M^e'^'^ one obtains for M^i 

This program has been realized for r decays and the extraction of a<j(M^) and ms{M^) 
in a series of papers [531 EH IS5] • 
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4.1 Resummation using the four- loop (5 function 

The resummation of the PT series for observables related to r decays with the four-loop 
MS /3 function has been studied in Ref. [63]. The integrals over the contour give the 
resummed functions that are analytic at the origin with a finite radius of convergence. 
The convergence radius of a series given by an integration over the contour and can be 
analyzed through the singularity structure in the complex plane [6^. In the lowest 
order example, the radius of convergence is determined by the solution of the equation 
1 — iTiPoa = 0, leading to a region of convergence given by \a\ < I/tt/^q. In higher orders of 
the P function the convergence properties of the resummed functions Mi^niO'r, become 
quite involved. The evolution of a along the contour = m^e**^, G [— tt, it] is governed 
by the renormalization group equation 

da 

- i— = f3(a) = -a^(l + ca + c^a^ + Csa^ + . . . ), (31) 

and the closest singularity in the complex plane then determines the convergence radius 
of the resummed functions Mj ^^(a^, The results for critical values of = ^^(m^) for 
increasing orders of the /3 function read 

= 0.444, af^ = 0.331, af) = 0.310, a^^^ = 0.299. (32) 

The convergence radii become smaller as the order of the /3 function increases. It is 
interesting to speculate about the possibility that the convergence radius continues to 
shrink as one goes to ever higher orders of the (3 function, including the possibility that 
the convergence radius shrinks to zero when the order of the perturbative /3 function 
expansion goes to infinity. 

The value of is outside the convergence region. This means that the perturbative 
approximation for the moments diverges at the scale determined by the experimental 
data for the semileptonic r-decay width. The resummed values are not accessible by 
using higher and higher order approximations of PT in polynomial form. 

4.2 Relations between observables using resummation 

Here we discuss a model of how the technique of using the direct PT relations between 
observables can give results that are also obtained in a more sophisticated resummation 
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approach |6l]. Consider two observables given by the perturbative series in some given 
scheme, i.e. 

f{a) = a{l-a + a^-...) = Y^ (33) 

and 

g{a) = a{l-2a + 4a^-...) = (34) 
The functions /(a) and g{a) can be seen to be related by 

9{f) = Y^ = f{l-f + f -■■■)■ (35) 

If we fit the right-hand side of Eq. (I33l) to an experimental value of about / = 0.6, we get 
a = 1.5. But for this value of the coupling, the series in Eq. ( I33l) diverges. One cannot 
get a from the series without an appropriate resummation procedure which, in this case, 
is trivially given by the exact formula. As a consequence we cannot get a prediction for g 
using the series in Eq. ( IMI) in terms of a. On the other hand, the direct relation in terms 
of the series in Eq. (135|) converges perfectly and gives an unambiguous result for g in 
terms of the measured /. Of course, for such an improvement to occur one has to analyze 
in detail the underlying theory and the origin of the series. The analysis of the r system 
taking into account corrections has been performed in finite order PT in Refs. [661167] . 
and with resummation in Ref. [68]. Resummation on the contour along the above lines 
in an effective scheme has been done in Ref. [65]. It happens that numerical results of 
the different techniques are different. It is important to understand the difference, or to 
answer the question to what extent the resummation recipe restores the same correlation 
function. 



5 Comparison of resummation techniques 

Because of the arbitrariness of resummation, it is important to understand the relation 
between different resummation techniques. Clearly, different resummations of asymptotic 
series give functions that have the same asymptotic expansion but differ from the original 
functions. General lore is that the difference between different resummation schemes 
behaves as exp(— l/a^). This form emerges from the Borel resummation recipe and is 
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confirmed by an explicit example with resummation on the cut [60] , although other forms 
are also possible [69] . 

Here we discuss the relation between CIPT and the resummation on the cut known 
as analytic PT [TQl [71] , following the lines of Ref . [72] . At leading order the moments can 
be expanded in a convergent series in for /3o«t < 1- Within the contour technique of 
resummation the finite radius of convergence is a general feature which persists in higher 
orders of the /3 function [HU |63] . The convergence radius decreases when higher orders 
of the /3 function are included. In practice, for af^ jii = 0.14 the relation a^^'P < l//3o is 
still approximately valid. The exact expression provides an analytic continuation beyond 
the convergence radius even when lies outside the convergence radius. 

We consider the moment moo and proceed with the analysis by constructing an efficient 
computational scheme. By integrating n times by parts, one obtains 



"^00 




n-l 



^ ^'^ vrr ; 2 \ TT J (1 + z/3oa,<^/7r)" ( ^ ' 



where the polar coordinate functions r and are defined by 



1 lb if3oar = re^""^, r = \J 1 + Plal, <t> = arctan(/3ottr)- (37) 

The n-fold integration by parts has removed a polynomial of order n from the expansion 
of the logarithm. 

One gets an asymptotic expansion where the residual term, i.e. the last term in 
Eq. f p6|) . is of the formal order a"^. However, the result is not a series expansion in 
the original coupling determined in the Euclidean domain but an expansion over a more 
complicated system of functions emerging in analytic PT [7T]. The system of functions 
is ordered and the asymptotic expansion is valid in the sense of Poincare. The system of 
functions is obtained by using the expression for the running coupling in the Euclidean 
domain and continuing it into the complex plane and onto the cut. When the analytic 
structure of the initial function is known, asymptotic expansions which converge fast for 
the first few terms (as a representation in the form of Eq. (l36l) ) are useful for practical 
calculations. The expansion in Eq. f p6|) can give a better accuracy (for some n and ar) 
than a direct expansion in Og. Indeed, this expansion includes a partial resummation 
of the TT^ terms which is a consequence of the analytic continuation [T7]. Therefore, the 
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expansion can be understood as being done in terms of quantities defined on the cut. 
Because the region near the real axis is important, the continuation causes a change of 



the effective expansion parameter — > a^/ + . The first term in the expansion 
shown in Eq. (!36|) is just the value for the spectral density expressed through the coupling 
in the Euclidean domain. 

With a concise expression for the moments at hand one can change the form of the 
residual term. The relation 

n — 1 ! 



TT 



;i + il3oar(p/Tr) 



is valid for any n and leads to a representation of the zeroth order moment of the form 

n-l 



moo = —r \ 



-/^°-+0 + X:(j-l)!(^)^in(j0) 
j=i V / 

n - l)\ f poarY ( r\ e'^dip 



J \J^oo Jn y (1 +i/3ottrV5/vr) 

Note that a "nonperturbative" term e~'^/*"^ has appeared. 



(39) 



The moments are analytic functions of a,- for small values of the coupling a^. This 
means that the non-analytic piece in Eq. (138|) cancels the corresponding piece in the 
residual term. If the residual term is dropped, the analytic structure drastically changes 
depending on which representation, either Eq. (136|) or ( 139|) . is used. 

One can recover the integral form of the moments as integrals over a spectral density 
by going to the complex plane in cp (see Fig. [1]) 

moo = 7r^ / TT~^ r ^0 

27r^ 1 + %l3oar(p/ir 

This representation is different from Eq. fl36l) for n = 1. The difference is an integral 
which can be explicitly computed, 

ttr dip 1 / n \ / 

TT~^ — = — arctan(/3oa^). (41) 

Now we consider the integration over a rectangular contour in the complex ip plane. 
The part of the contour on the real axis from — vr to it leads to the moment expressions. 
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Figure 1: Integration contour in the complex (p plane 



The integral over the contour is given by the residue at the pole ip = in/Poar- We thus 
have 

With the substitutions -7r//3o«r = \n{A^/M^), = \n{s/M^) one obtains 



Finally 



1 /A2\ 1 /-^^^ / ln(s/A2) , _ 

^00 = TTo + / arccos , " ' tTo • (44) 




7r2 + ln"(s/A 

One recognizes this representation as an integration over the singularities of n(g^). In 
addition to a cut along the positive semi-axis, there appears also a part of the singularity 
on the negative real s axis. This part is a pure mathematical feature of the definite 
approximation chosen for 11 (g^). The result reads 

with 

a{s) = ^6{A^ + s)9{-s) + -^6{s) arccos | _H^/£^^ \ ^^g^ 
^0 VyWh?(^)/ 

This formal result can be rewritten as an integration over the spectrum a{s) using 
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Cauchy's theorem. Indeed, 

Disc U{s) = '^le{A^ + s)e{-s) + -e{s) arccos | _^(^/£)^ | y (47) 

which coincides with (t(s) in Eq. fj46l) . The part of the spectrum on the positive real axis 
is an analytic continuation of the function n(Q^) to the cut [T71 UHl |18l |70]. It can be 
written in the form 

'^c(-5) = ^arctan(/3oa(s)). (48) 
ttPo 

The differential equation determining the continuum part ads) through its initial value 
cr(M^) can be constructed by differentiating Eq. PHI) with respect to s, 

^<M£)_ /oWy 1^ (49) 

as \ n J 1 + /3Qa(s)^ 

By inverting Eq. fHHj) one has /Soal-s) = tan(7r/3oO"c(s)) for s > 0. Therefore, one obtains 

s-^ads) = ^sin^(7r/3o(Tc(s)) for s > 0. (50) 

as Po 

This equation can indeed be considered as an evolution equation for the spectral density 
(Tc{s) determining a^s) through its initial value ct(M^). Therefore, one can introduce an 
effective charge aM{s) = nads) with a corresponding evolution equation 

s — -. = T^sm (7r/3oaM(s)). (51) 

as po 

Thus, one defines the coupling as the value of the spectral density on the cut far from 
the infrared region. The evolution of this coupling is calculated by taking into account 
analytic continuation. In this case, the coupling has an infrared fixed point with the value 
flA/(0) = l//3o- If Adler's function starts with another power of the coupling constant, as 
is the case for gluonic observables, this picture will change. For 

D(Q')^(^^y (52) 

the spectral density in the leading-order /3-function approximation reads 

[il ln^(s/A2) + TT^ 7r2(l + Plads)) ' 
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Figure 2: Contours in the complex plane starting and ending in the region A. The different 
contours talce into account possible occurrences of singularities. The left hand side of the 
figure shows the standard circular path which circumvents the singular region B. Other 
singular regions as discussed in the text (regions B' and B") lead to different possibilities 
for choosing a path (Ca, Cb, Cc). The path Cd crosses the singular region and, therefore, 
cannot be used from the perturbation theory point of view. 



and an effective coupling is 



<^m{s) = " (54) 



The /3 function for the effective coupling obtained from Eq. (j5 

Puiau) = -/^oOmv^I - (vr/3oaM)^ = -/^oOm ^1 - ^(tt/^oOm)^ + ^(a^)^ (55) 
differs from the (3 function of the effective coupling oa/ (c.f. Eq. (15T]) ). 

fiuiaM) = ^ sin^ (vr/3oaAf) = -f^oalj ( 1 - -^(vr/JoaA/)^ + 0{alj)] (56) 

at next-to-leading order. Thus, resummation on the contour and on the positive semi-axis 
of the s plane differs by the integral over the negative real semi-axis for s. 

In the context of the contour formulation it is not essential what particular point-by- 
point behaviour exists in the infrared region. For the analytically continued correlator this 
is not important unless the contour crosses a nonanalytic region. Whatever singularities 
exist in the infrared region (Regions B, B', or B" in Fig.[2]), the contour includes them. The 
resummation on the contour is explicitly perturbative. For the resummation on the cut. 
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the extrapolation of the running of the couphng constant to the infrared region is crucial 
since one has to interpret the integration over the infrared region. Formal manipulations 
with t = Poa{Ml) \vl{M'^/s) give 

"^'^ ~Jo 1 + /3o«(M2) ln(./M2) - /3o X l-t ' ^ ^ 

Not all expressions in ( 1571) are well-defined. This is particularly true for the third form 
being a Borel representation. The problem can be reformulated as a divergence of the 
asymptotic series. Indeed, by expanding the expression for the running coupling under 
the integration sign in a PT series one has 

"''„W,,.^„,(^)". (58) 

The summation of the series in Eq. f l58|) is related to the interpretation of the integral. 
Therefore, an integrable behaviour of the coupling constant at small s offers a recipe 
for the summation of the asymptotic series. This solution is strongly model dependent 
because the extrapolation of the evolution into the infrared region is essentially arbitrary. 
The explicit form of the extrapolation in Eq. (146|) gives an extrapolation motivated by 
analytic continuation. It can also be considered as a special change of the renormalization 
scheme [60] . Indeed, for the coupling qm with the evolution given in Eq. (15T1) one obtains 
an infrared fixed point. Without the negative part of the spectrum that emerges in the 
exact treatment of quantities originally defined in the Euclidean domain the analytic 
moments can be defined by 

•""" = /*"^ (59) 

One can relate these moments to the moments on the contour by the relation 

^or' = -^e--/*"^+moo (60) 
Po 

which contains an explicit "non-perturbative" term. 

The question of uniqueness of the resummation is important. However, it is only 
relevant if many terms of the PT expansion are known. It happens that this uniqueness can 
be usefully employed in multi-loop computations [731 IZl]- Note that it is also important 
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to have examples of higher order PT behaviour in simple QFT models as unique record- 
breaking computations show [75l [76]. In some cases uniqueness is good laboratory for 
checking multi-loop techniques [771 EHl [79]. The methods are also applied in some exotic 
nonlinear sigma model, topological theories, and SUSY [80] . 

6 Summary 

We have reviewed different ways of interpreting perturbation theory results for the de- 
scription of r-decay observables. Experimental data are very precise and theory matches 
it by unprecedented numbers of PT terms. Higher orders of PT are available for many 
cases that allows both for the description of data and for the extraction of the parame- 
ters of the theory with high accuracy. However, the PT series converge slowly requiring 
improvement that can be achieved through: 

i) manipulation with schemes for a set of related observables avoiding artificial inter- 
mediate quantities such as MS quantities; 

ii) different resummation techniques; 

iii) the avoidance of unnecessary expansions - treating final polynomials in the coupling 
as exact expressions, e.g. /3 function and anomalous dimensions in renormalization 
group equations for coefficient functions of operator product expansions 

and some other improvements that you name. This richness is available due to RG 
properties that allows one to control the scaling behaviour and the invariance of the 
theory. 
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A Appendix 

The phenomenology of hadronic r decays is contained in the correlator of the weak cur- 
rents jj^{x) = cos(6'c)m7m(1 - 75)c? + sm{ec)u'j^{l - 75)5, 



with p{s) = Im n'^^^(s + i0)/7r and 
The total r-decay rate 



p{s)ds 
s — q^ 



is a useful observable which is measured with high experimental precision. Since p^^'^{s] 
is a distribution one considers moments of the spectral density of the form 



_ (fc + / + l)! /-^^^Y . V ( s \'p{s)ds _ 
km Jo V J \M? / 



which contain information about the hadronic spectral density p^^'^{s). The function p{s) 
is related to Adler's function by 

p{s)ds 



where = —q^ and D{Q'^) are computable in PT. In the mass zero hmit, the PT 
expression for Adler's function reads 

D{Q'^) = l + as + ha] + fego^ + fcga^ + O(a^) 

with as = as{Q'^)/'n- The MS-scheme coefficients read 

299 , , 58057 779,, , 75,, , , 

fc, = 3^-9C(3), I, = _ - _C(3) + -C(5) (Al) 
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while /cs 



= 49.08 [391 SI]. The correction to the total width 5p is 

5^^ = a, + 5.2023a2 + 26.366a^ + (78.003 + 49.08)a^ + O(a^) 



which corresponds to the experimental value 5'^^^ = 0.216 ± 0.005. 
The renormalization group equation for a(/_i^) reads 

^ ^ -a2(/3o + Pia + /^aa^ + Psa^ + . . . ) 

with 

^ n . n 3863 ^ 140599 445,, , 
A = -, ft=4, = — , = — + -C(3) 

for = = 3 [8I1I19]. 

The coefficients /j„ are given by 

fon = /(0,n), /i„ = /3o/(l,n), = /3o' (/(2, n) + pi/(l, n)) , 

/sn = /3o'(/(3,n) + (/,(2)- 1,(1)2 -i7r2)/(l,n) + ^Pi/(2,n)+p2/(l,^) 

/4. = /3o'(/(4,n)-3(/,(2)-/,(l)2-l7r2)/(2,n) 
+2(J.(3) - 3/,(l)/.(2) + 2/^(1)3) J(l, n) 

13 ~ 1 ~ _ _ 

+Pi(-/(3, n) + 5(/,(2) - 1,(1)2 _ -vr2)/(l, n)) + 3p2/(2, n) + p3/(l, 



with 



J(m,n) = . , , /,(m) = 2J(m,0)-2J(m,2) + J(m,3), 
(n + 1)^ 



m 



/(m,n) = I{m) + ^y\lr{v)I{ 

p=0 



m — p,n). 



The Pi are scheme-independent quantities given by 

Pi = p2 = -L[l32- I3iki + (3o{k2- kl)], 

Po Po 

P3 = ^[/33-2f32k^ + f3,kl + 2f3o{k3-3kik2 + 2kl 
Po 
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